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Abstract 

Pair production processes of neutral Higgs particles will allow us to study the 
trilinear Higgs couplings at future high-energy colliders. Several mechanisms 
give rise to multi-Higgs final states in hadron interactions. In the present 
paper we investigate Higgs pair production in gluon-gluon collisions. After 
recapitulating pair production in the Standard Model, the analysis of the cross 
sections is carried out in detail for the neutral Higgs particles in the minimal 
supersymmetric extension. 
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1 Introduction 



The reconstruction of the Higgs potential is an experimental prima facie task to establish 
the Higgs mechanism as the basic mechanism for generating the masses of the fundamental 
particles. This task requires the measurement of the self-couplings of the Higgs particles. 

In the Standard Model (SM) [0] the trilinear and the quartic couplings of the physical 
Higgs particle H are uniquely fixed if the Higgs mass is known. 

In the minimal supersymmetric extension of the Standard Model (MSSM), a large 
variety of couplings exists among the members h, H, A, H ± of the Higgs quintet 0. [h 
and H are the light and heavy CP-even Higgs bosonsQ, A is the CP-odd (pseudoscalar) 
Higgs boson, and H ± is the charged Higgs pair.] While general two-doublet models 
contain three mass parameters and seven real self-couplings in CP conserving theories, 
the Higgs self-couplings are fixed in terms of gauge couplings in the MSSM, and the mass 
parameters can be expressed by the two vacuum expectation values of the neutral Higgs 
fields and one of the physical Higgs masses. Since the sum of the squares of the vacuum 
expectation values is given by the W mass, tg/3, the ratio of the two vacuum expectation 
values, and Ma, the mass of the CP-odd Higgs boson A, are generally chosen as the free 
parameters of the MSSM. The trilinear and quartic Higgs self-couplings are determined 
by those two parameters. 

The measurement of the Higgs self-couplings will be a very difficult task. In the 
Standard Model the production cross section for HH Higgs pairs are small, similarly the 
continuum production of Higgs pairs in the MSSM. Only if heavy MSSM Higgs bosons 
can decay into pairs of light Higgs bosons, the associated Higgs self-couplings can be 
determined fairly easily by measuring the decay branching ratios. Several aspects of 
multi-Higgs production have been discussed in the literature. Most analyses treat HH 
pair production in the Standard Model only at the theoretical level 0, |[]. Within the 
MSSM, the search for the heavy Higgs boson H has been simulated in the bb decay mode 
at the LHC ||. h — > AA events have been searched at LEP ||; more general aspects of 
multi-Higgs production in e + e~ collisions have only recently been analyzed theoretically 
in Ref.@. 

Several mechanisms give rise to the production of pairs of neutral Higgs bosons in 
hadron collisions. Multi-Higgs final states can be produced through Higgs-strahlung off 
W bosons and through WW fusion in proton-proton collisions. Bremsstrahlung of Higgs 
particles off heavy quarks can also be exploited. In the present paper we discuss the 
production of Higgs pairs in gluon-gluon collisions at the LHC. We have determined the 
cross sections for the continuum in the Standard Model 

VV-> 99 -> HH (1) 

as well as for the continuum and resonance decays in the minimal supersymmetric theory 

pp^gg^ ^i<t>j = h,H,A (2) 

1 Following standard notations we have taken care that no confusion will arise from using the same 
symbol for the SM and the heavy CP-even MSSM Higgs particle. 
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restricting ourselves to neutral Higgs bosons in the present report^]. It can be anticipated 
from the large number of low- x gluons in high-energy proton beams that the gg channel 
is of particular interest in the continuum for fairly low Higgs masses. In the MSSM the 
on-shell production of heavy Higgs bosons with subsequent Higgs cascade decays will 
eventually provide a copious source of light multi-Higgs final states. The cross sections 
are affected, besides the normal Higgs-boson couplings to gauge bosons and fermions [[/J , 
by the trilinear Higgs couplings M: 

SM : \huh = ^§- (3) 



MSSM : X hhh = 3 cos(2a) sin(/3 + a) 



3e cos 3 a 



Af| sin/3 

^Hhh = 2 sin(2a) sin(/3 + a) — cos(2a) cos(/3 + a) 



3e sin a cos 2 a 



Ml sin (3 

» _ . , , n \ s . , n x 3e sin 2 a cos a 
*HHh — — 2 sm(2a) cos(p + a) — cos(2a) sm{p + a) + — j 



^hhh = 3 cos(2a) cos(/3 + a) + 



Mf sin (3 
3e sin 3 a 



Mf sin/3 
e cos a cos 2 (3 
Mj sin/3 

e sin a cos 2 /? 
Mf sin/3 



A/iaa = cos(2/3) sin(/3 + a) + 

\haa = - cos(2/?) cos(/3 + a) + — ^ — ^ (4) 



The couplings in the SM as well as the MSSM are normalized to Ao = \\I~2G p^l 2 M\. The 
MSSM couplings depend on /3 and the mixing angle a 

Ml+MZ 

tan 2a = f tan 2/3 (5) 

M\ -M| + e/ cos 2(3 H K ' 

in the CP-even Higgs sector; radiative corrections have been included in the leading m\ 
one-loop approximation, parametrized by 



3Gp mj 
V27T 2 sm z p 



Ml 

! + -f 

mi 



(6) 



with the common squark mass fixed to M$= 1 TeV. In our numerical analysis we have 
included the leading two-loop corrections to the MSSM Higgs masses and couplings, taken 
from Ref.||. The behaviour of the couplings with Ma is shown for two representative 
values of tg/3 = 1.5 and 30 in Fig.l. 



2 The derivation of the SM cross sections has nicely been described in Ref.Q. We will briefly discuss 
this case again to exemplify the techniques for the simplest scenario of Higgs pair production, before 
generalizing the analysis for unequal masses and pseudoscalar Higgs particles in the MSSM. 
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The paper is organized as follows. In the next section HH pair production through 
gluon collisions will be discussed in the Standard Model. In the third section we will 
present the results for all pair combinations of neutral Higgs bosons in the MSSM. Phe- 
nomenological aspects including Higgs-strahlung and WW fusion will be presented in a 
sequel to this paper [KJ. 



2 Higgs Pair Production in the Standard Model 

Two mechanisms contribute to the production of Higgs pairs through gg collisions in 
the Standard Model, exemplified by the generic diagrams in Fig.2a/b. (i) Virtual Higgs 
bosons which subsequently decay into HH final states, are coupled to gluons by the usual 
heavy-quark triangle [|ll|, (ii) The coupling is also mediated by heavy-quark box 
diagrams. 

In the triangle diagram Fig. 2a the gluons are coupled to the total spin S z = along 
the collision axis. The transition matrix element associated with this mechanism can 
therefore be expressed by the product of one (gauge invariant) form factor Fa, depending 
on the scaling variable tq = imq/s, and the generalized coupling Ca defined as 

r - a M| (7) 

C A - X HHH j-jjf (7) 

The coefficient Xhhh denotes the trilinear self-coupling HHH in the Standard Model, 
cf. eq.(§). s is the square of the invariant energy flow through the virtual Higgs line. The 



well-known form factor [JIT, [12] is given in Appendix Al. 

The box diagrams in Fig. 2b allow for S z = and 2 gluon-gluon couplings so that the 
transition matrix element can be expressed in terms of two (gauge invariant) form factors 
Fa and Gn- Fairly compact expressions of the form factors are given in Appendix Al 
where the Higgs masses have to be specified to M c = = Mh- The couplings between 
the Higgs bosons and the quarks are normalized to unity by definition, 

C D = 1 (8) 

[The overall normalization is included in the prefactors of the cross section and the form 
factors.] The form factors obtained in this way for the simplified case of equal masses 
agree with the results in Ref. M. 

In the two limits of light and heavy Higgs bosons with respect to the loop quark mass, 
very simple expressions can be derived for the three form factors Fa and Fa, Gn- 

(i) Large Quark Mass Limit 

The form factors can be evaluated either by taking the limit tuq ^> s ~ M\ in 
the Feynman amplitude or, equivalently, by exploiting the elegant low-energy theorem 
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Fa = ™>Qd(F&/rnQ)/drnQ [12-^] for external light scalar Higgs bosons. This reduces the 
complexity of the calculation considerably: 



Fa = 3 + O (s/ml 



and 



Ga 



-\ + 0(s/ml) 
O (s/m 2 Q ) 



(9) 



(10) 



2 form factor 



Both S z = form factors F& and F n survive in this limit, whereas the S z 
G n vanishes asymptotically. 

(ii) Small Quark Mass Limit 

In the limit of light quark masses compared with the invariant energy^, large logarithms 
of the light quark mass log s/ttiq will occur in the form factors. Detailed inspection of 
the Feynman amplitudes leads to the following final expressions: 



rn 



Q 



log ( J + m 



+ °\f~ 



and 



Ga 



O (ml/s 
O (ml/s) 



'IV 



(12) 



All form factors vanish to 0(mn) in this limit since they are suppressed by the Yukawa 
coupling oc nig; this is well-known for the quark triangle. [This limit is not of much 
practical use in the SM; yet it will be relevant later in the MSSM where 6-quark loops 
are dominant for large tg/3.1 



The differential parton cross section can finally be written in the form 
da(gg - HH) 



G 2 F a 2 s 



dt 



256(2tt) 



Ca Fa + CaF n \ + \CaGa\ 



(13) 



where t is the momentum transfer squared from one of the gluons in the initial state 
to one of the Higgs bosons in the final state. The total cross section for HH Higgs 



3 In the leading logarithmic approximation the scale is set effectively by the Higgs mass since the 
invariant energy y/I is of the same order as Mh for the cross sections we shall study. 
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pair production through gg in pp collisions^ can be derived by integrating (|1^) over the 
scattering angle and the gg luminosity 

a{pp -v gg-> HH) = / dr——a(s = rs) (14) 

jam 2 h /s ar 

This cross section has been evaluated numerically. The analysis has been carried out for 
the LHC cm. energy y/s = 14 TeV; the top quark mass has been set to m t = 175 GeV 
|[L5|| . The result is shown in Fig. 3. For SM Higgs masses close to Mjj ~ 100 GeV the cross 
section is of order 10 fb. However, it drops rapidly with increasing Higgs mass; this is a 
consequence of the fast fall-off of the gg luminosity with rising r [which is of the order 
of x 2 g in the pp collisions]. Nevertheless, at Mh ~ 100 GeV order 2,000 events will be 
produced for an integrated luminosity of / £ = 100/6 , giving rise to 46 and bbrr final 
states with large transverse momenta. 



3 Neutral Higgs Pairs of the MSSM 

A large variety of neutral Higgs pairs of the MSSM can be produced in gluon-gluon 
collisions: 

PP^ 99 -> hh, hH, HH, hA, HA, AA (15) 

The analysis of the cross sections is in general much more involved than in the SM for 
two reasons. First, the masses of the two Higgs bosons in the final state are in general 
different, and second, besides pairs of CP-even Higgs bosons, also pairs of CP-odd and 
mixed CP-even/CP-odd pairs can be produced. 

a) Pairs of CP— even Higgs bosons hh, hH, HH 

The mechanisms for the production of CP-even MSSM Higgs bosons gg — ► HiHj 
(Hi = h, H) are similar to the Standard Model. The generic triangle and box diagrams 
are shown in Figs.4a/b. All possible combinations of light and heavy Higgs bosons h and 
H can be coupled in 3-particle vertices. 

The matrix element of the triangle contributions can be expressed in terms of the form 
factor Fa, given in Appendix Al, and the couplings 

C A = Ci + C% (16) 

with 

cr^»^ m s _ KiH f iMhiHThiH 9r (in 

The couplings <?q denote the Higgs-quark couplings in units of the SM Yukawa coupling 
[v^2G , i?] 1//2 mQ. They are collected in Table [I]. The Higgs self couplings XiiiHAh/H) can be 
read off eq.(f|). 

4 In analogy to single Higgs production via gluon fusion, QCD corrections are expected to enhance the 
event rate significantly. However, they are known only for the triangle and not for the box; therefore 
they are not taken into account in the numerical analysis. 
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9t 9b 


SM H 


1 1 


MSSM h 
H 
A 


cos a/ sin (5 — sin a/ cos /3 
sin a/ sin /3 cos a / cos /? 
l/tg/3 tg/3 



Table 1: ffiggs couplings in the MSSM to fermions relative to SM couplings. 



For the box contributions we obtain analogous expressions, with two independent form 
factors F a and G u contributing. Their expressions can be found in Appendix Al. The 
corresponding couplings C a are are built up by the Higgs-quark couplings, 

Cu = g$9Q i (18) 

The two limiting cases in which the Higgs masses are either much smaller or much 
larger than the loop quark masses, are both physically relevant in this case since the top 
quark as well as the bottom quark play a role, with M| -C mf and M\ 3> vn^. The limits 
of the form factors are the same as for the SM Higgs bosons eqs.(|9HT2|). 

From the differential parton cross section, which includes t and b quark loops, 

~ ] ~d HiHi) = J|^[ |E(Ca Fa +^^)I 2 + IE^Go| 2 ] (19) 



t,b t,b 



the pp cross section can be derived by integrating over the scattering angle and attaching 
the gluon luminosity The result of the numerical analysis is shown for the final states 
hh, hH and HH in Figs.5a/b/c for two representative values tg(3 = 1.5 and 30. [Mixing 
in the stop sector is not taken into account.] 

The final state hh deserves special attention. The cross section for moderate values 
Mh is large since the heavy CP-even Higgs boson H can decay into the two light Higgs 
bosons (dashed lines). H is produced directly in gg collisions, coupled through the quark 
triangle, so that the matrix element involves one power of the electroweak coupling less 
than continuum Higgs pair production. Except for a dip due to a zero in the Hhh 
coupling Xnhh, the cross section in the resonance range is of order 10 3 to 10 4 /6, dropping 
to the typical hh continuum value of 30 fb for increasing H Higgs masses. For tg/3 =1.5 
the curve shows a sharp threshold behavior for Mh ~ 2m t due to on-shell top pair 
production in the dominant top quark triangle and, moreover, due to the sharp fall-off 
of the branching ratio BR(H — > hh). The cross section for tg/3 = 30 depends strongly 
on the Higgs boson masses, because resonance production gg — > H — ► hh is kinematically 
forbidden for 110 GeV < < 210 GeV. Above this range this channel opens up again 
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with a small branching ratio ~ 10 3 giving rise to the small bump in the cross section at 
M A ~ 210 GeV. 

The cross sections for hH and HH final states are much smaller, cf. Figs.5b/c. This 
is the result of the absence of any resonance effect and the phase space suppression due to 
the large Higgs masses in the final states. The cross section for pp — > gg — > hH depends 
strongly on tg/3 and it varies between very small values and ~ 100 fb if the entire Higgs 
mass range is sweeped. For tg/3 = 1.5 the Higgs mass dependence is smooth, whereas for 
tg/3 = 30 a sharp peak develops at M A ~ 100 GeV due to the strong variation of the 
MSSM couplings in this region. At large Higgs masses the cross section decreases strongly 
with increasing tg/3. A similar picture emerges for the pair production of heavy scalar 
Higgs bosons pp — > gg — > HH. Its size increases at large Higgs masses for increasing 
tg/3, where the bottom quark loops become dominant. The strong variation of the cross 
section for M A ~ 100 GeV and tg/3 = 30 is again due to the rapidly changing MSSM 
couplings in this mass range. 

b) Mixed CP-even/CP-odd pairs 

The analysis of the CP-mixed hA and HA final states is in many aspects different 
from the previous cases. First of all, the final state can be produced through the decay of 
a virtual Z boson as shown in the generic diagrams of Fig.6. By evaluating the Feynman 
diagrams one finds the following generalized charges [Hi = h, H] 



(i) s-channel A exchange: 



M 2 

Ci = W * s - M\ + iM A T A $ (20) 



(ii) s-channel Z exchange: 



M 2 

° l = XzAH > § - M\ +%M Z T Z aQ (21) 

where a<g = 1(— 1) for top (bottom) quarks denotes the axial charge of the heavy 

loop quark Q [note that only the Z axial coupling of the heavy loop quark can 
contribute], and the trilinear couplings XzAHi are given by 

XzAh = -cos(/3-a) (22) 
^zah = sin(/3 - a) 



(iii) box contributions: 

= 9 A Q 9 H q (23) 

The expressions of the form factors F^j Z and F n , G a can be found in Appendix A2. 

In analogy to the CP-even case, simple expressions can be derived for the form factors 
in the limits of large and small Higgs masses compared with the quark masses. 
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(i) Large Quark Mass Limit 



The triangular form factor can be derived in this limit from the axial anomaly. The box 



form factor is given by the derivative of the anomaly F a 
Simple expressions can be obtained: 



and 



A 



Gn 



1 + O (s/m 2 Q 
O (s/ml 



-1 + [s/m l Q 
0(s/m 2 r 



m 2 Q d(F£/m Q )/dm Q [0,0. 



(24) 



(25) 



(ii) Small Quark Mass Limit 



F n 
Gn 



m 



log I — — I + in 



n2 + o(^ N 



s - Ml M* t - M\ 



M\ 



1 + 



m, 



log 



m, 



Q 



O {m%/ 8 
O (ml/s 



+ m 




(26) 



The differential parton cross section is determined by the form factors, including t and 
b quark loops: 

da<S3 7 t AHi> =^H^m C - F - + CoFo)|* + |£Co G o|=] (27) 



t,b t,b 



The results of the numerical analysis of the two pp cross sections are shown in Figs.7a/b. 
For light masses M^ the invariant mass [hA] of the final state is small so that the cross 
section is enhanced by the large phase space. For large Higgs masses the cross section 
decreases strongly with increasing tg/3 and drops to a level of about 10 _1 to 10 3 /6 at 
Ma ~ 1 TeV, where it cannot be observed anymore. The strong variation for tg/3 = 30 
at Ma ~ 100 GeV is generated by the rapid change of the MSSM couplings with M4. 

Due to the larger mass of the heavy scalar Higgs boson H, the cross section for HA 
production turns out to be smaller than for hA production. 
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c) AA pairs of CP— odd Higgs bosons 



This case is again closely related to the CP-even pairs. In the box diagrams, Fig. 8b, 
some terms flip just the sign when the 75 x 75 couplings are reduced to II in the loop trace. 
As a result, the matrix elements can be expressed in terms of the couplings 



C A = Ci + Cf (28) 



with 



ih/H , Mz h/H 



G A = Xaa(h/h) 



M 2 /H + iM h/H T h/H 



Cu = (29) 

The form factors F& for the triangle contributions and F a , Gu for the box contribution 
can be found in Appendix A3 by choosing M c = Md = Ma- 

In the limits of large and small Higgs masses compared with the quark masses simple 
expressions can be derived. 

(%) Large Quark Mass Limit 



The triangle form factor coincides with the scalar expression eq.(|]), and the box 
form factor can be obtained from the derivative of the gluon self-energy [Ai(ggA 2 ) = 
\f2GpmQdM. [gg)/dmQ\ |14]] leading to the results 



A 



1 0(s/m%) (30) 



3 



and 



Fa = l + 0(s/m 2 Q ) (31) 
G Q = 0(s/m 2 Q 



(ii) Small Quark Mass Limit 



m. 



log ^ 



27T 



o 



m 



Q 



and 



F u = 0(m 2 Q /s) 
G n = 0(m 2 Q /s 



(32) 



(33) 
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The pp cross section is shown in Fig. 9. The large value of the cross section for small 
A masses is due to resonance H — > AA decays. When this channel is closed, the cross 
section becomes more and more suppressed with rising A mass in the final state. For 
tg/3 =1.5 the cross section develops a kink at the (tt) threshold Ma = 2m t due to S- wave 
(if)-resonance production. Above this mass value the phase space suppression leads to 
a rapidly decreasing cross section, to a level of ~ 10~ 3 fb at Ma ~ 1 TeV. For tg/3 = 30 
the bottom quark loops are dominant so that the cross section depends smoothly on Ma- 
For large pseudoscalar masses the signal increases with increasing tg/3. 



4 Appendices 



Al. Form factors for two scalar Higgs bosons g a gb — > H c H d 

notation: 



parameter definitions: 

8 = (Pa+Pb) 2 , 

S = s/m 2 



t = {Pc ~ Pa) 2 , U= (p c - Pbf 

T = t/m 2 , U = u/m 2 



p c = M 2 c /m 2 Q} p d = M 2 /m 2 Q} r Q = A/S 



Ti = T - p c , 
Scalar integrals: 



u x = u- p c , 



T 2 =T - p d , U 2 = U - p d 



27T 2 ( g 2_ m 2:) ( g+p .)2_ m | (q+ Pi + Pj ) 2 -ml 



D 



ijk 



d 4 q 



1 



(q 2 - m 2 Q ) (q + p t ) 2 - m 2 Q (q + p, L + pj) 2 - m 2 Q (q + p { + pj + p k ) 2 - m 2 Q 



The analytic expressions can be found in Ref . [ 16 

triangle form factor: 



Fa = ^{2 + (A-S)m 2 Q C ab } = r Q 



1 + (1 - T Q )f{T Q ) 



f(r Q ) 



arcsm 



log 



T Q 



1 2 



TQ>1 
TQ < 1 
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box form factors: 



F n = ^^S + 8Sm 2 Q C ab -2S(S + p c + p d -8)m 4 Q (D abc + D bac + D acb ) 

TiC ac + UiC bc + U 2 C ad + T 2 C bd - (TU - p c pd)m 2 Q D acb J 
(T 2 + p cPd - 8T)m 2 Q \sC ab + T x C ac + T 2 C bd - STm 2 Q D bac 



G 



+ {pc + Pd 

1 



S(TU - pcpd) 



+ (U 2 + PcPd ~ 8U)m 2 Q SC ab + U x C bc + U 2 C ad - SUm 2 Q D abc 

_( T 2 + c/2 _ 2pcPd )(T + [/ - 8)mJC cd 

-2(T + [/ - 8)(TC/ - p cPd )m 4 Q (D abc + D bac + D acb ) j 



tensor basis: 



S? = 



S z = 2 



PaPb 
(PaPb) 



A »v _ n V , P 2 cPaP b 2(PbPc)P&Z 2(PaPc)P b Pc ^P 



with 



PliPaPb) PriPaPb) 

pI 



PriPaPb) Pi 



2 _ (PaPc)(PbPc) 2 

p T — Z, 



(PaPb) 

A 1 -A 2 = 0, A 1 ■ A 1 = A 2 ■ A 2 = 2 



matrix element: 



M(g a g b ^H c Hd) = M h A + M% + M n 



mT 



Gf^sS h/H p . u v r 

2^/2 n A ^a^i^ e a e fe d ab 
Gpct s s 



C a (F n A lflu + G a A 2 ) e^e v b 5, 



>ab 



A2. Form factors for a mixed scalar-pseudoscalar pair g a gb — > A c H d 

triangle form factor: 



K = 



-2m\C ab = T Q f(r Q ) 

S \ Pc~ Pd 



M 2 Z/ 



S 



1 + 2m 2 Q C ab 



Ml 



Pc - Pd 

S 



1 - T Q f(T Q ) 
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box form factors: 



F a = -gjj -2S{S + p c - p d )m i Q (D abc + D bac + D acb ) 



+ (pc-pd)m 2 Q 
1 

S(TU - p c p d ) 



T\C ac + U x C hc + C/ 2 C ad + r 2 C M - (TU - p c/ o d )mjD ac6 ] J 



(U 2 - p c pd)m\ 



SC ab + U x C hc + £/ 2 c7 ad - SC/mjD, 



aba 



-(T 2 - p c p d )ml 



+ 



'S'Cab + T'lCac + T 2 C b d — STrr?qD bac 
(T - EOmJC* 



+2(T - C/)(TC/ - p cPd )m%(D abc + L> boc + D ac6 )| 



tensor basis: 
S z = 



A? 



AS? 



c P-Vp a p b 



(PaPb) 

pH e VPaPbPc _|_ e fJ.paPbPc _|_ Lp^p\ e Hup a Pc _|_ (j) a p c )e IJ ' UPbPc 



(PaPb)pj 



with 



(PaPc)(PbPc 



Pi 



(PaPb) 

A x ■ A 2 = Ai • Ai = v4 2 • A 2 = 2 



matrix element: 



A^i + .Mf + .Mo 

2^2^ A A ^ a 6 a 
Gpa s s 

2^ 



>ab 



A3. Form factors for two pseudoscalar Higgs bosonsf] g a gb — > A c A d 

triangle form factor: 

2 



5 



{2 +(4 -5)771^}=^ 



1 



;i - r Q )f(r Q ) 



5 We present the formulae for different masses of the two pseudoscalar particles so that they can also 
be applied to non-minimal SUSY extensions in which several pseudoscalar particles in general belong to 
the Higgs spectrum. 
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box form factors: 



F n = ^^AS + 8Sm 2 Q C ab -2S(T + U)m 4 Q (D abc + D bac + D acb ) 

+ (p c + pd)m 2 Q TxCac + UiC bc + U 2 C ad + T 2 C hd - (TU - p c pd)m 2 Q D acb J 
(T 2 + p c p d )m 2 Q SC ab + TiC ac + T 2 C bd - STm 2 Q D bac 



G 



1 



S(TU - pcpd) 



+ (U 2 + p c p d )m^ SC afe + U x C bc + C/ 2 C ad - SUm 2 Q D abc 
_( T 2 + f/2 _ 2pcP(i )(T + U)m 2 Q C cd 

-2(T + C/)(TC/ - p cPd )m%{D abc + D bac + L> ac6 )} 



tensor basis: 
S z = 
S z = 2 



PaPb 
(PaPb) 



AP, , PlPaPb 2(PbPc)lfi% 2(PaPc)P b Pc , 2p£p, 



with 



PriPaPb) PriPaPb) PriPaPb) Pt 

9 _ n (PaPc)(PbPc) 2 

Pt — z 7Z~ZT\ Pc 



(PaPb) 

A 1 ■ A 2 = A x ■ A 1 = A 2 ■ A 2 = 2 



matrix element: 



M(g a g b ^A c A d ) = M h A + Ml + M n 
U h/H _ G F a s s h/H 



ab 



Gpa s s 
2^2ix 



C n (FnAi^ + GuA 2jiU ) 



>ab 
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FIGURE CAPTIONS 



Fig.l: Trilinear Higgs couplings in the MSSM as a function of the pseudoscalar Higgs 
mass Ma for two representative values tg/3 =1.5 and 30 [no mixing in the stop sector, 
c.f. Ref.[|]. 

Fig. 2: Generic diagrams describing SM Higgs pair production in gluon-gluon collisions: 
a) triangle and b) box contributions. 



Fig. 3: Total cross section of SM Higgs pair production at the LHC [y/s = 14 TeV]. The 
top mass has been chosen to be m t = 175 GeV and the bottom mass = 5 GeV. The 
GRV parametrizations of the parton densities have been adopted. The factoriza- 



tion and renormalization scales are identified with the invariant mass of the Higgs boson 
pair. The world average value a s (Mz) = 0.118 [13 has been chosen for the QCD coupling. 



Fig. 4: Generic diagrams contributing to pair production of CP-even MSSM Higgs bosons 
in gluon-gluon collisions gg — > hh, hH, HH: a) triangle and b) box contributions. 



Fig. 5: Total cross sections for pair production of CP-even MSSM Higgs bosons in gluon- 
gluon collisions at the LHC: a) hh, b) hH and c) HH pair production for tg(3 = 1.5 
and 30. The secondary axes present the scalar Higgs masses M h / H corresponding to the 
pseudoscalar masses Ma for both values of tg/3. The same parameters and parton densi- 
ties as in Fig.3 have been adopted. The dashed curves in Fig. 5a represent the resonant 
contributions gg — > H — > hh. 

Fig. 6: Diagram of the Z boson s-channel contribution to mixed CP-even/CP-odd MSSM 
Higgs pair production in gluon-gluon collisions and box diagram. 



Fig. 7: Total cross sections for mixed CP-even/CP-odd pair production of MSSM Higgs 
bosons in gluon-gluon collisions at the LHC for tg/3 = 1.5 and 30: a) hA and b) HA 
production. The secondary axes show the scalar Higgs boson masses corresponding to 
the pseudoscalar mass Ma for both values of tg/5. The same parameters as in Fig.3 have 
been chosen. 



Fig. 8: Generic triangle and box diagrams for pair production of CP-odd MSSM Higgs 
bosons in gluon-gluon collisions. 

Fig. 9: Total cross section of CP-odd MSSM Higgs pair production in gluon-gluon col- 
lisions at the LHC for tg/3 =1.5 and 30. The parameters are the same as in Fig.3. The 
dashed curves represent the resonant contributions gg — > H — > A A. 
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Fig. 4 




Fig. 5a 
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Fig. 5c 
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Fig. 6 




Fig. 7a 
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Fig. 7b 
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